Atomic multiplet calculation of 8^5/2 4/ resonant x-ray diffraction from Ho metal 
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We compare for Ho metal the x-ray absorption spectrum and the resonant soft x-ray diffraction 
spectra obtained at the 3ci5/2 — > 4/ {K'h) resonance for the magnetic 1st and 2nd order diffraction 
peaks (0, 0, r) and (0, 0, 2r) with the result of an atomic multiplet calculation. We find a good 
agreement between experiment and simulation giving evidence that this kind of simulation is well 
suited to quantitatively analyze resonant soft x-ray diffraction data from correlated electron systems. 

PACS numbers: 75.25.-(-z, 78.20.-e, 78.70.Ck, 78.70.Dm 



Resonant soft x-ray diffraction (RSXD) at the transi- 
tion metal L2,3 (2p — > Sd) or lanthanide M4^5 {id 4/) 
resonance is a relative new tool to study orbital, charge 
or magnetic ordering phenomena^^'^'^'^'^'^'^'^ '^^'^^i^^'^'^ 
The technique makes use of the high sensitivity of an elec- 
tronic transition at one of these resonances to the elec- 
tronic and magnetic state of the scattering ions, which 
renders RSXD directly sensitive to super structures in- 
volving spatial variations of the electronic state. 

In many interesting cases like transition-metal L2.3 
RSXD from transition-metal oxides and A'h^^ RSXD from 
lanthanide systems the resonance involves a strong inter- 
action between the core-hole and the valence electrons 
making the resonant diffraction spectra strongly exci- 
tonic. A simple picture of RSXD in terms of the density 
of unoccupied states is hence not appropriate. In order 
to understand RSXD data from these systems quantita- 
tively one needs a more realistic theoretical description. 
Fortunately, such a description has been developed for 
the closely related technique of x-ray-absorption spec- 
troscopy (XAS) at transition metal ^2,3 edges and rare- 
earths M4,5 edges Jiii^ii^iii As RSXD and XAS involve 
the same optical transitions, i.e., the same orbitals are 
involved and the same transition matrix applies, one can 
expect that the theoretical treatment used for XAS can 
directly be apphed to RSXD. 

For the description of the excitonic intermediate state 
in RSXD, i.e., the final state in XAS, one needs to include 
the full electron-electron interaction. For transition- 
metal oxides this is usually done within a cluster cal- 
culations including the transition-metal ion and its oxy- 
gen ligands.^^'^^'^^ In lanthanide systems the overlap be- 
tween the 4/ wave function and its surrounding is usu- 
ally so small^* that an atomic model (possibly includ- 
ing a crystal field) is sufficient. For XAS this technique 
has been used to obtain invaluable information on the 
local electronic structure of many transition metal and 
rare earth compounds. For RSXD, on the other hand, 
there are up to now only a few calculated spectra avail- 
able in the literaturei^ i^'^'^^i^° Although these calcula- 



tions match experimental data fairly well, they describe 
rather complex systems such that material properties and 
the applicability of the methods itself are probed at the 
same time. In order to verify that a description in anal- 
ogy to XAS data is really appropriate for RSXD spectra 
a test with a simple model system is desirable. 

In this paper we present a comparison of theoretical 
and experimental RSXD and XAS data from Ho metal 
at the M5 threshold. Ho has a particularly simple and 
well established magnetic structure and has served as a 
reference system for the development of resonant mag- 
the conventional^^'^^ and soft 
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x-ray 

Ho crystallizes in the hep-structure; we refer 
in the following to the unit cell spanned by three orthog- 
onal lattice vectors with a — 3.58 A, 6 = 6.20 A, and c = 
5.62 A. Below Tn = 131.2 K bulk Ho orders with a heli- 
cal magnetic structure (Fig. [T]) . The propagation vector 
of the helix r lies in the c direction. The spin direction is 
perpendicular to Resonant scattering from this heli- 
cal structure reveals two Bragg peaks (Fig. [2), with one 
of them at (0, 0, r) corresponding to the total periodicity 
of the helix and the second one at (0, 0, 2t) corresponding 
to half its periodicity.^'* In the resonant magnetic scat- 
tering amplitude derived by Hannon et al.^ the (0, 0, r) 
peak is described by the term that is linear in the direc- 
tion of the magnetic moment, m, while the (0, 0, 2r) peak 
relates to the term that is quadratic in m. 

At the Ho M5 resonance the photon penetration depth 
changes strongly with energy. The corresponding change 
of the scattering volume leads to distortions of the reso- 
nance spectra when obtained from bulk samples^. For 
the experiment we therefore chose a 11-monolayer (ML) 
thin Ho film for which the probed volume is defined by 
the film thickness rather than by the photon-penetration 
depth. In the sample the Ho layer is sandwiched be- 
tween Y layers with the consequence that the helix pe- 
riod length in this film of about 7.5 ML^^ is shorter than 
the bulk value of about 10 ML. The resulting larger r 
value in the film allows for a better separation between 
magnetic peaks and specular reflectivity background at 



FIG. 1: Schematic representation of the helical magnetic 
structure of holmium metal. The large arrows within the 
grey box symbolize the magnetization direction within the ab 
planes. 



low momentum transfer, q. T^v is reduced with respect 
to the bulk value to about 76 K.^^ While bulk Ho shows 
at low temperatures a lock-in transition to a commensu- 
rate magnetic structure and a canting of the spins out 
the ab plane to a conical phase we observed none of 
these effects in our film. 

The diffraction experiments were carried out at the 
soft x-ray beamline UE46-PGM1 at the electron storage 
ring BESSY using the soft x-ray diffractometer designed 
at the FU Berlin. The experimental RSXD data were 
obtained by recording the scattered intensity at the re- 
spective peak maximum as a function of photon energy, 
keeping the momentum transfer constant. XAS data 
were taken from Ref. and were obtained from a 31 
ML thick Ho film grown in situ on a W(llO) substrate. 
The cluster calculations have been done with the program 
XTLS8.3r^ developed by A. Tanaka and previously used 
for the description of RSXD on Lai.8Sro.2Ni04ii The lo- 
cal model used for these calculations is the same as ap- 
plied to describe the XAS on a series of Lanthanides by 
Thole et al.^ It includes the full electron-electron re- 
pulsion, whereby we scaled the 4/ — 4/ Slater integrals 
to 65% and the 3d — 4/ Slater integrals to 75% of the 
Hartree-Fock values. The Hartree-Fock values have been 
obtained with the use of Cowans programi^ The calcu- 
lations have been done in spherical symmetry and on a 
single configuration. This is not a necessary restriction, 
as extensions to crystal and ligand field theory can easily 
be madeJ^ii^ii^iii For Ho metal the crystal-fields acting 
on the 4/ shell are small, though, and therefore spherical 
symmetry is sufficient. 

XAS and RSXD can be most easily formalized if one 
starts from linear response theory. For each atom one 
can define a scattering tensor Fj which describes the re- 



FIG. 2: L-scan through the (0, 0, r) and (0, 0, 2r) magnetic 
peaks at 26 K. The increasing intensity at low L-values is 
caused by specular reflectivity from the sample. 

lation between the components of the scattered radial 
wave at a distance R of the scatterer and those of 
the incoming plane wave E^^ at frequency to. El^f{uj) — 
Fj(uj)-El^{uj){e''^'^l''/R). The scattering tensor at a spe- 
cific scattering vector q can be written as a sum over the 
scattering tensors of each atom multiplied by the appro- 
priate phase: F{q,uo) = e*''"'^" where 

the first sum runs over all unit cells at positions i?„ in 
the probed volume and the second sum over the ions at 
positions Tj in one unit cell. If one now defines Sin (Sout) 
as the polarization of the incoming (outgoing) photons 
one can write the intensities of RSXD and XAS as: 

lRSXD{q,i^) = \el^f ■ F{q,uj) ■ e„i\'^ (1) 
Aire 

IXAS (^) = 3 [4. • F{q = 0, C^) • £,„] (2) 

The relation between the absorption and scattered inten- 
sity is given by the optical theorem, which follows from 
conservation of energyi^ The real and imaginary part of 
F{q, u) are related by the Kramers-Kronig transforma- 
tions in the uj domain. 

The scattered intensity depends on the polarization of 
the incoming and outgoing photons; hence it is a 3 x 3 ten- 
sor, whereas the absorption depends only on the incom- 
ing polarization and is described by the diagonal tensor 
elements only. By rotating the polarization vector (real 
as well as complex rotations) one can write all compo- 
nents of the F tensor as a linear combination of diagonal 
terms. For example, defining (x+y) as the light polarized 
in the x + y direction and {x + ly) as right circular polar- 
ized light, we can describe the scattering of s-polarized 
incoming to y-polarizcd outgoing light by: 

1 + « 

Fxy — -F'(a:+y)(2;+i/)+«-F'(2:+ij,)(j:+ij,) - {Fxx+Fyy) (3) 

This gives a direct relation between XAS spectra 
recorded with different incoming light polarizations and 
all possible RSXD channels. 
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FIG. 3: (color online) Experimental and theoretical Ho M5 
XAS and RSXD spectra. 



In Fig. [3] we show the theoretical and experimental 
XAS and RSXD spectra of the Ho M5 edge. Although 
the probed resonance is the same for all three experi- 
ments, one can clearly see that each measurement, XAS, 
RSXD at g = T, and RSXD at g = 2r has a very dif- 
ferent spectral shape. The theoretical calculations repro- 
duce these different spectra quite well. The agreement 
between theory and RSXD is equally good as the agree- 
ment between theory and XAS. This is pleasing as it is 
known that cluster calculations can reproduce XAS for 
many transition metal and rare-earth compounds rather 
well. Obviously the level of agreement is the same for 
our RSXD spectra. 

The huge differences in spectral shape between XAS 
and RSXD at different q values can be understood if one 
realizes which different elements of the F((f, oj) tensor are 
probed by the different signals. We chose a notation with 
respect to an orthogonal coordinate system aligned with 
the crystalline axes. For a magnetic ion with the magne- 
tization in the x direction in otherwise spherical symme- 
try Fj{uj) can be written as: 



^ XX 


TP 

^yx 


F,x \ 














^yy 


F.y 







F± 


iFq 


Fxz 


Fy. 


F.. J 




^ 


-iFq 


F± 



(4) 



The diagonal elements of Fj{uj) refer to scattering pro- 
cesses which do not alter the polarization of the scat- 
tered light. The off-diagonal elements describe scatter- 
ing processes involving changes of the light polarization. 
Magnetic scattering breaks time-reversal symmetry and 
is hence (on a basis with one of the principle axes parallel 
to the magnetization, in the present case x) described by 
an antisymmetric off-diagonal part of Fj{ijj). 

It is tempting to relate the symmetric part of the F 
tensor to charge scattering and anti-symmetric part to 



magnetic scattering. Care should be taken, however, as 
a local magnetization also produces symmetric contribu- 
tions to the F tensor. One reason is that a local magneti- 
zation can induce an orbital momentum, which is related 
to a non-spherical charge distribution and therefore gives 
rise to charge scattering. Birt there is also an effect creat- 
ing symmetric contributions to the F tensor due to a local 
magnetization that is not related to a local charge dis- 
tribution, but pure magnetically: At L2.3 or M^^^ edges 
the core hole is created in a shell with orbital momentum 
(2p or 3c?) and large spin-orbit coupling. The magnetiza- 
tion of the valence electrons interacts via the multipolar 
part of the Coulomb interaction (direct and exchange) 
with the core electrons. This results in a shift of res- 
onance energy between different polarizations, reflected 
in symmetric and antisymmetric contributions to the F 
tensor. This second contribution would for example also 
be present in or configurations. 

In order to analyze the different contributions to the 
different spectra in Fig. 1 we use the simplest sufficient 
model, which is reproducing the two observed Bragg- 
peaks. It consists of a chain of four atoms labeled A,B,C, 
and D stacked in the z direction with the magnetic mo- 
ment alternating in the x, y, —x, —y direction. In our 
model we have four different scattering tensors. Fb is 
related to Fa by a 90° rotation around z and so on. 
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The generalization to helices with longer periods or in- 
commensurate order is straightforward. The neglect of 
crystal field effects and assumption of a spherical sym- 
metry are no necessary simplifications, but for Ho turns 
out to be sufficient, most probably because crystal field 
effects are smaller than spin-orbit coupling and magnetic 
interactions. 

The F tensor corresponding to a phase repetition after 
four lattice spacings [q — (0, 0, r)] then becomes: 

F[q = (0, 0, t),u] = Fa + iFb ~Fc~ iFd 




-2Fq -i2Fq 




(6) 



The F tensor corresponding to a phase repetition after 
two lattice spacings [q = (0, 0, 2t)] is: 

F[q - (0, 0, 2t),uj] =Fa-Fb + Fc-Fd 
2{F\\ -Fi_) 

2(F^-F||) I (7) 











The (0, 0, r) peak is hence only given by the off-diagonal 
element of the single-scatterer's tensor, Fq, and the 
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FIG. 4; Experimental geometry. The direction of the Poyint- 
ing vectors are determined by the Bragg condition. The polar- 
izations can be defined as: a — (0, 1,0), nin = (sin 0, 0, cos 0), 
and TTout = (— sin 8, 0, cos 6). 

(0, 0, 2r) peak only by the difTerences between the di- 
agonal terms F± — i^|| . This property is also conserved 
for other helix periodicities. 

The F tensor for XAS [q= (0,0,0)] becomes: 

F[q ^ {0,0,0),uj]^Fa + Fb + Fc+Fd 

/ 2(^^11+ Fx) \ 

2(Fi + F||) (8) 

\ 4Fl y 

In our case the XAS measurement have been preformed 
in the paramagnetic phase above the Neel temperature 
which means that the local spin-direction should be av- 
eraged over all possible directions and the F tensor for 
this case becomes: 

F[q = (0,0,0],c^) = 

/ F|| + 2Fi \ 

= 4/3 F|| + 2Fj_ (9) 

\ F|| + 2F_L / 

We have so far discussed the scattering tensors within 
the reference frame of the crystal, which makes the math- 
ematical formalism a bit more transparent. Usually one 
writes the polarization vectors with respect to the direc- 
tions parallel to the scattering plane (tt) and perpendic- 
ular to it (cr) . For Ho with the scattering vector in the z 
direction and if we choose a to be in the scattering plane 
we find the polarizations to be equal to: 

<J^n^'Jout = (0,1,0) (10) 

TTm — (sin6', 0, cos6') (11) 
TTout = (-sin6', 0,cos6') (12) 



with = oj/2 and cu « 25° for (0,0, r) and w « 51.4° 
for (0, 0, 2r) in the considered energy range. With the 
scattering angles as defined in Fig. 4. One now has all 
ingredients to straightforwardly calculate the RSXD and 
XAS intensity. By substituting the F tensor as given in 
formula 6,7, and 9 into equation 1 or 2 for RSXD and 
XAS respectively one finds the intensity with the polar- 
izations as given in formulas 10-12. Note that since the 
F tensors as written in equation 6, 7, and 9 depend only 
on a single function, the spectral lineshape does not de- 
pend on the polarization or on the rotation of the sample 
around q. The total intensity still does, though. 

The analyzes of the scattering intensity and line-shape 
as presented here can be directly related to the results 
derived by Hannon et al^, if one changes the basis of 
the F tensor from linear polarized light {x,y,z) to circular 
polarized light (-1 = l/\/2(x-iy),0 = z,l = l/V2{-x- 
ly)). For a site with a magnetization in the z direction 
one finds that 1/2{F( + FiJ = 1/2(F4 + F^y) = F^, 

l/2{Fl, - Fi) = l/2z(F4 - FJJ = F^, and Fo = F,,. 
The formulas of Hannon et at for dipole interactions 
are retrieved if one uses this basis and realizes that in 
spherical symmetry the F tensor for a magnetization in 
arbitrary direction is given by a simple rotation of the 
F tensor magnetized in the z direction; Fj(6, (j)) ~ Re.4, ■ 
Fj{z) • i?g 0. Whereby Fj{z) is the local F tensor for 
a magnetization in the z direction, Fj [9, (p) the local F 
tensor for a magnetization in the 9, direction and 
a rotation matrix rotating a vector in the z direction to 
a vector in the 0, <j) direction. 

To summarize, we have shown that it is straightfor- 
wardly possible to simulate RSXD spectra from a 4/ sys- 
tem using configuration interaction cluster calculations 
with very good agreement between experiment and simu- 
lation. This technique is well developed for XAS and can 
without further modification directly be used for RSXD 
spectra. With the use of such analyzes one can expect to 
obtain a wealth of spectroscopic information from RSXD 
on the ordered electronic structure like orbital, spin and 
charge ordering quantitatively. 
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